Chiral excitonic order from twofold van Hove singularities in kagome metals

Recent experiments on kagome metals AV3Sb5 (A=K,Rb,Cs) identify twofold van Hove singularities (TvHS) with opposite concavity near the Fermi energy, generating two approximately hexagonal Fermi surfaces – one electron-like and the other hole-like. Here we propose that a TvHS generates a novel time-reversal symmetry breaking excitonic order – arising due to bound pairs of electrons and holes located at opposite concavity van Hove singularities. We introduce a minimal model for the TvHS and investigate interaction induced many-body instabilities via the perturbative renormalisation group technique and a free energy analysis. Specialising to parameters appropriate for the kagome metals AV3Sb5, we construct a phase diagram comprising chiral excitons, charge density wave and a region of coexistence. We propose this as an explanation of a diverse range of experimental observations in AV3Sb5. Notably, the chiral excitonic state gives rise to a quantum anomalous Hall conductance, providing an appealing interpretation of the observed anomalous Hall effect in kagome metals. Possible alternative realisations of the TvHS mechanism in bilayer materials are also discussed. We suggest that TvHS open up interesting possibilities for correlated phases, enriching the set of competing ground states to include excitonic order.

Kagome systems have been a major focus of theoretical and experimental investigation; due to their ability to realise Dirac points, flat bands and van Hove singularities, they have been predicted to host a range of novel correlated phases of matter [1][2][3][4][5][6][7][8][9] . Recently, a new class of materials AV 3 Sb 5 (A=K,Rb,Cs) have attracted a great deal of attention due to their demonstration of unconventional superconductivity alongside competing density wave order, spatially modulated superconducting order and possible signatures of Majorana states in superconducting vortices  . Unusually, the materials exhibit timereversal symmetry breaking with an anomalous Hall conductivity in spite of the absence of magnetic ordering; the origins and relationship between superconductivity, competing order, and the anomalous Hall effect remain an open question.
The materials consist of a stack of two dimensional layers-a kagome lattice of vanadium and antimony alternating with a hexagonal lattice of antimony and triangular lattice of the alkali metal K/Rb/Cswith electrical transport predominantly in-plane, as demonstrated by the large ratio between the out-of-and in-plane resistivity R c /R ab ≈ 600. The Fermi surface of these materials consists of several distinct contours, including nearly circular contours centred at the Γ and K points as well as two approximately hexagonal contours 44 . Systems with hexagonal Fermi surfaces, corresponding to saddle-points in the electronic dispersion, have been predicted to give rise to chiral superconductivity and competing density wave order, due to the effects of Fermi surface nesting 45 . However, ARPES and DFT results reveal that the hexagonal Fermi surfaces in the vanadium metals exhibit an unusual feature-twofold van Hove singularities (TvHS), for which the saddle-points at each Fermi surface possess opposite concavity, resulting in one electron-like Fermi surface and one hole-like Fermi surface 39 .
We argue that doping a system to a TvHS has an ineluctable influence on the low-energy physics. A single vHS results in a tendency towards density wave ordering and superconductivity. The appearance of TvHS introduces an additional tendency towards excitonic order-corresponding to a condensation of electron-hole pairs-owing to the coupling between an electron-like and hole-like Fermi surface. We introduce a low-energy model, which incorporates the TvHS-featuring an electron-like and hole-like Fermi surface, each doped near their respective vHS. To understand the interplay and competition between the various many-body instabilities, we employ the perturbative renormalisation group (RG) method to determine the dominant ground state order [45][46][47][48][49] , complemented by a Landau-Ginzburg free energy analysis of competing ground states. A chiral excitonic order naturally emerges, which breaks time-reversal symmetry and exhibits a quantum anomalous Hall effect. The chiral excitonic state appears as a generic weak coupling instability, but explicit modelling for AV 3 Sb 5 suggests these materials exist in an intermediate coupling regime; guided by ab initio results we generate a phase diagram featuring charge density wave order, chiral excitonic order and a region of coexistence. We suggest that the phenomenology encompassed by the TvHS model accounts for key features observed in the vanadiumbased kagome metals, and could further motivate TvHS engineering in van der Waals heterostructures and bilayer materials.

Tight-binding Hamiltonians with TvHS
A TvHS consists of two Fermi surfaces with opposite concavity vHS so that one surface is electron-like and the other hole-like-e.g., arising from doping near the M-point of a 2D hexagonal Brillouin zone as shown in Fig 1a. The opposing concavities of the respective saddlepoints can be seen from the colour plot in Figs. 1c; going from outside to inside the hexagonal Fermi surface, the energy changes sign, but the sign change is opposite for the two Fermi surfaces. The two Fermi surfaces may arise due to two hexagonal (honeycomb or kagome) bilayers, or a single layer with two sets of orbitals-the latter case is the origin of the TvHS in vanadium metals AV 3 Sb 5 .
To be explicit, we will introduce a particular lattice model realisation of a TvHS. A tight-binding model of a kagome monolayer with two sets of orbitals that has been used to describe AV 3 Sb 5 is given by where a y i,ν creates fermions on site i and in orbital ν = c, d. The differing orbital potentials, ϵ c − ϵ d ≈ t c + t d shift the energies of the two bands, aligning their valence and conduction bands and resulting in a TvHS.
The bandstructure of a TvHS can be realised in both honeycomb and kagome systems (we discuss alternative tight-binding models in the Supplementary Material). However, the orbital content of the wavefunction at the M-points is qualitatively different in these two cases. For honeycomb, with two sublattices, the wavefunction at the Mpoints has equal support on both sublattices. Meanwhile, in kagome systems the wavefunction near the M-points exhibits the so-called 'sublattice interference effect' 8 : at a given M-point, the conduction band wavefunctions have support only on one sublattice and are referred to as p-type (owing to their 'pure' sublattice composition) while the valence band wavefunctions have support on the other two sublattices and are referred to as m-type (due to their 'mixed' sublattice composition). The sublattice structure has important consequences when considering interaction effects, as we discuss below.

Patch model
The problem of interaction driven instabilities on a single hexagonal Fermi surface (i.e., single vHS) has been previously studied using a three patch model 45 , whereby the full Brillouin zone is restricted to three momentum space patches near the vHS at the M-points, since they dominate the density of states. Following this approach, we define a three patch model and further introduce a flavour degree of freedom to account for the two opposite concavity Fermi surfaces-fermions of one flavour (created by c † ) are electron-like, while fermions of the other flavour (created by d † ) are hole-like, The patch index α = 1, 2, 3 indicates a fermionic excitation within a cutoff Λ of the momentum M α . Setting μ = 0 corresponds to doping exactly to the TvHS. The opposite concavity is encoded by the relative minus sign between the c and d dispersions. The TvHS should be contrasted with the problem of fixed concavity vHS with multiple flavours 50 -the opposite concavity of the two vHS plays a fundamental role in the interaction physics. The patch dispersion take the saddle where t ν is a characteristic energy scale, and equals the nearest neighbour hopping of the ν-fermions (ν = c, d) in the simple tight-binding model. Fermions at patches α ≠ β are connected by the nesting vector Q αβ = M α − M β , for which ε ν k + Q αβ ,β ≈ À ε ν k,α . Making contact with ab initio results for AV 3 Sb 5 , the c-(p-type) and d-(m-type) fermions arise from the vanadium d yz and d xz orbitals respectively, and have t c ≈ 0.5 eV, t d ≈ 1 eV 51 . In the patch model, this sets t c /t d ≡ κ = 2. For completeness we will analyse both κ = 1 and κ = 2. It is known from ARPES that the c-band vHS is near-perfectly nested, while the d-band vHS exhibits quartic corrections 39 . Close to the Mpoint these corrections are subdominant to the quadratic part of the dispersion, and hence only influence the ultraviolet behaviour of the theory, near the cutoff Λ ≈ 0.5 eV. Since our analysis probes infrared scales far below Λ, it is well-justified to ignore the quartic corrections.
Below we will analyse three distinct cases: (i) honeycomb systems, for which the sublattice support on the two-flavour vHS is the same, in the particle-hole symmetric limit κ = 1; (ii) kagome systems, in which the two-flavour vHS have different sublattice support, i.e., m-and ptype, with κ = 1; (iii) kagome systems with κ = 2, which we have argued to describe kagome metals AV 3 Sb 5 .

Interactions
We now consider the possible couplings between the fermions. Owing to the large density of states near the TvHS the Coulomb repulsion is expected to be strongly screened and we therefore model the interactions as short-ranged. The most general set of interactions between patches/flavours allowed by momentum conservation are where V g,ν are intraflavour couplings, V h are interflavour densitydensity couplings, V j are flavour pair hopping, and V l are flavour exchange couplings, resulting in 20 independent interactions. A schematic illustration of the g, h, j, l couplings, as well as their representation in terms of Feynman diagrams, is shown in Fig. 1c, d. Additional details are found in the Supplementary Material. In the kagome case, projecting the sublattice wavefunctions onto the Coulomb interaction results in different intraflavour couplings depending on whether the flavour has pure or mixed sublattice structure, a manifestation of the sublattice interference effect in kagome patch models. We have therefore allowed for different couplings V g,ν on each flavour. Performing this projection explicitly and using the calculations of ref. 51 gives the estimates of the bare coupling values shown in Table 1. The values taken from ref. 51 are defined at the lattice scale; using these as input to our effective theory neglects the renormalisation flow between the lattice scale and Λ. The sublattice interference effect has crucial consequences for the bare couplings: for instance, on a p-type vHS, the wavefunctions at different patches are orthogonal, and so the interpatch Coulomb repulsion is suppressed, resulting in g 1,c = g 3,c = 0. Thus, any attractive interactions present in the system, for e.g., due to phonons, immediately result in attractive couplings.
In the honeycomb case, the orbital form factors are the same for both flavours and so we expect V g,c ≈ V g,d . This reduces the number of independent coupling constants from 20 to 16. We shall present results for both models below.

Instabilities
Considering the interacting Hamiltonian, we determine which instabilities arise within the framework of RG. The instability of the metallic phase and onset of an ordered ground state is signalled by the susceptibility of the associated order parameter: the strongest ordering tendencies are those with most divergent susceptibility. In the case of a nested Fermi surface, a density wave instability arises because the nesting condition ε p ≈ − ε p+Q , implies the total energy of a particle with momentum p and hole with momentum p + Q is approximately zero. Similarly, the energy of an electron and a hole at the TvHS is ε c p + ε d p ≈ 0. Without including interactions, it costs zero energy to create either of these particle-hole states, and hence, for an arbitrarily small attraction between particles and holes the system becomes unstable to lowering its energy by spontaneously creating many such pairs, analogous to the usual superconducting instability. The RG method provides an unbiased approach to study competing orders on an equal footing, by resumming the logarithmically divergent corrections to the bare couplings and determining which ordering tendency dominates [45][46][47][48][49] .
In Table 2, we enumerate the ordered states which naturally arise in the TvHS model, i.e., those with nesting tendencies. The first three-CDW, SDW and SC-occur in the case of a single vHS. The next threesinglet and triplet excitonic order, as well as interflavour pair density wave (PDW)-are new instabilities introduced by the TvHS.

RG analysis
We turn now to the RG treatment which identifies the leading instabilities, i.e., the dominant ground states in Table 2. Firstly, we compute the leading log 2 corrections to the bare couplings defined in (3). The equations define how the couplings evolve with the RG time t, which is a proxy for the energy scale; here t → ∞ corresponds to taking T → 0. The full RG equations for our model are lengthy, since they involve twenty independent interaction constants (Fig. 1d), so we state their general form here and reserve explicit expressions for the Supplementary Material. The RG equations describing the flow of the couplings g i , h i , j i , l i (where i = 1, 2, 3, 4) take the form Projecting the pure and mixed sublattice form factors onto the cRPA results of ref. 51  where β g i,ν , β h i , β j i , β l i are functions of all twenty couplings. Secondly, we compute the leading log 2 corrections to the order parameters, which generates the linear set of gap equations, Diagonalising the gap equation matrix V ij and integrating over the RG time t, one identifies the leading eigenvalue λ i (t), which diverges fastest with t. The associated eigenvector is the order parameter with the largest critical temperature T c = Λe À1=ðν 0 λ i Þ 1=2 , and is therefore the dominant order at T ≲ T c . Multiple orders of comparable T c may arise, in which case one must compute the Landau-Ginzburg free energy to ascertain whether such phases compete or coexist.

Dominant instabilities
A subset of the couplings diverge with increasing RG time t → ∞. In this limit, the diverging couplings tend towards fixed constant ratios of each other referred to as fixed rays. The relative magnitudes of the couplings determine which ground state dominates. All possible choices of bare initial coupling values flow to one of these possible sets of ratios in the deep infrared, which therefore represent universal properties of the model. We now present the set of fixed rays possible in our TvHS patch model (for a derivation see the Supplementary Material). Despite the large number of interaction terms there turn out to be only a small set of fixed rays, shown in Fig. 2 Crucially, in all cases d-wave excitons emerge at a fixed trajectory, demonstrating the naturalness of excitonic order. As in the case of single vHS 45 , we find that d-wave superconductivity is also a natural instability of the TvHS model.
For arbitrarily small initial couplings, the fixed rays are reached at long RG times, which corresponds to the deep infrared. However, the initial couplings could be sufficiently large that an instability occurs before the fixed ray is reached. In such a case it is appropriate to instead explicitly compute the flow from a specific set of initial conditions, and examine when an instability is reached. Given the significant magnitude of the bare values of the couplings in AV 3 Sb 5 (Table 1), we believe that such an analysis is more appropriate when comparing with experiment, and is presented below.

Properties of the chiral excitonic condensate
We focus attention on some key properties of the excitonic phases which appear. The two d-wave excitonic order parameter structures which appear are Here, Φ a , Φ b are real scalars, ϕ a , ϕ b are distinct U(1) phases, and the row vectors enumerate patch indices, therefore encoding the spatial structure of the two order parameters. Continuing to the full Fermi surface, the spatial vectors schematically behave as ∼ cosð2θ k Þ, sinð2θ k Þ, with θ k the momentum angle. Similar d-wave eigenvectors appear for the superconducting states Δ, which are the two-flavour analogues of the superconducting states found in 45 .
Near the critical temperature, the Landau-Ginzburg free energy is found to be where F 0 is the free energy of the free fermions, with expansion coefficients a Φ , c Φ > 0. The free energy is minimised by coexisting order parameters, with Φ a = Φ b = Φ 0 and ϕ a − ϕ b = π/2 (mod π). The coexisting states form a single order parameter of the form Φ = Φ 0 e ± iθ α , θ α = {a, b, c}. Continuing around the Fermi surface, the combined order parameter becomes Φ ∼ Φ 0,k ðcosð2θ k Þ ± i sinð2θ k ÞÞ, which is a chiral d ± id order. The chirality ± is spontaneously selected by the ground state, which therefore breaks time-reversal symmetry. In addition to broken TRS, the chiral order parameter winds twice along the Fermi surface and vanishes away from it, thereby exhibiting a non-trivial topology with Chern number |C| = 2. In order to illustrate this, we diagonalise a mean-field Hamiltonian (Methods III) defined on a lattice-we consider a two-orbital kagome lattice model in an infinite ribbon geometry with zigzag edges. The 1D dispersion of the ribbon is plotted in Fig. 3 for the d + id phase, which exhibits two chiral edge modes, with the left/right-movers propagating along the top/bottom of the ribbon. Full details of the lattice model are provided in the Supplementary Material. The non-trivial topological invariant implies a quantised anomalous Hall conductivity σ xy = Ce 2 /(2π), which is carried by two chiral edge modes. We note that this value of σ xy = Ce 2 / (2π) accurately accounts for the intrinsic contribution to the anomalous Hall effect seen in AV 3 Sb 5 17 .
Notation: α, β index patch, σν act on flavour, indexed by Latin characters ν = c, d with

Coexistence of chiral excitons and charge order
The RG procedure determines which phases are dominant, but is not sufficient to determine the actual ground state when two phases have comparable T c . Motivated by experiments on AV 3 Sb 5 , we examine the ground state when CDW and chiral excitons are proximate instabilities.
To carry out the analysis, we consider the Landau-Ginzburg free energy for chiral excitons and CDW, written as F = F Φ + F C + F Φ,C , with F Φ from Eq. (8) and The expansion coefficients a ν , b ν , c iν depend on temperature and are computed in the Supplementary Material. Here V νν 0 is the CDW gap equation matrix (6). Henceforth, we denote the largest eigenvalue of V νν 0 by λ CDW . In the free energy (9), there are six complex numbers, C αβν , describing CDW order. Physically the C αβν correspond to the magnitude of the order for the three distinct vectors Q αβ ∈ {Q 12 , Q 23 , Q 31 }, on the two distinct Fermi surface flavours. From the gap equation (6) we find that the leading CDW order has C αβc = ρC αβd , where ρ is a real number. In particular, |ρ| = 1 in the particle-hole symmetric limit of κ = 1. Moreover, the gap equation distinguishes real charge density order (rCDW) whereby C * αβν = C βαν and purely imaginary order (iCDW) whereby C * αβν = À C βαν . We separately considered parameter regimes in which rCDW and iCDW were the leading CDW phase.
We turn now to the phase diagram predicted by (9). To construct the phase diagrams of Fig. 4(a)i and 4(b)i we allow the eigenvalues λ CDW and λ Φ to be free variables. To illustrate the property of coexistence and our phenomenological proposal for these AV 3 Sb 5 , we take a realistic range of coupling eigenvalues, consistent with DFT calculations 51 , and fix T = 80 K (which enters via the coefficients a ν , b ν , c iν ). We find three distinct phases: (i) chiral excitons, (ii) CDW, and (iii) coexistence of chiral excitons and CDW. In the rCDW/iCDW phases, the 3Q state is favoured, i.e., CDW order is nonzero for all three nesting vectors Q αβ , and corresponds to C 12ν = C 23ν = C 31ν ≠ 0. In the excitonic phases, the chiral (TRS breaking) d + id state is favoured. In the region of coexistence, chiral excitons and the 3Q CDW are favoured. We point out that d-wave excitons coupled to a nematic CDW (e.g., C 12ν > C 23ν = C 31ν ) was observed as a local minima, but did not appear as the global minimum over the parameter range searched. Experiment indicates that TRS breaking and CDW coexist in AV 3 Sb 5 , and set in at T * ≈ 100 K. We propose that the coexistence phase demonstrated by our Landau-Ginzburg analysis provides a phenomenological explanation of the physics of kagome metals at T ≲ T * .

Truncated RG flow and phase diagram for AV 3 Sb 5
To complement the analysis leading to Fig. 4(a)i and 4(b)i, we now directly compute the eigenvalues λ CDW and λ Φ from the RG procedure. Unlike for the fixed ray analysis, here we must provide initial conditions for the RG flow. Once initialised, we perform the RG flow down from a UV scale of Λ ≈ 0.5 eV to an infrared scale set by T. We use the resulting renormalised couplings as input to the free energy, minimising to obtain the resulting ground state. This procedure generates the phase diagrams of Fig. 4(a)ii and (b)ii.
We discuss now the choice of initial couplings that lead to Fig. 4(a) ii and (b)ii. Given that several of the couplings in Table 1 vanish, we allow for the situation where these couplings take negative values. By inspection of the gap equation (6), we see that an initial value of g 1,i < 0 promotes CDW (this was first noted in 52 for the problem of a single vHS), while h 1 < 0 promotes chiral excitons. To this end, we first allow for both g 1,c , h 1 < 0, and subsequently arrive at the phase diagram of Fig. 4(a)ii. In addition, we allow for g 3,d = 0, and arrive at Fig. 4(b)ii. The phase diagram is qualitatively the same for g 3,d < 0. Next we mention that the magnitudes and ratios have been estimated from ab initio calculations (Table 1). A more accurate treatment would account for the renormalisation of the couplings in going from lattice to the patch UV cutoff Λ. We have not included these effects in our analysis. We stress that our use of the values in Table 1 is to illustrate that there exist physically reasonable bare couplings, which produce the desired phenomenology.

Discussion
We introduced and analysed a minimal model to describe interacting fermions near a twofold van Hove singularity (TvHS)-two opposite concavity vHS near the Fermi level. We found the opposite concavities of the two vHS crucially affect the possible many-body instabilities, relative to the single vHS case. In particular, excitonic order contends as a possible instability and generically results in a chiral d-wave excitonic phase in hexagonal systems such as honeycomb and kagome lattices. We contrast our scenario with topological excitonic states, which have been previously explored theoretically [53][54][55][56] ; in our case, the topology of the d + id condensate is not inherited from the Berry curvature at the K-points or from spin-orbit coupling, but appears at the M-point intrinsically due to interaction driven, spontaneous time-reversal symmetry breaking. These findings suggest a new class of candidate materials for topological excitonic ground states.
TvHS were recently seen experimentally in AV 3 Sb 5 39 . We now discuss key features of experiment and the extent to which the TvHS minimal model explains them: First, signatures of time-reversal symmetry breaking, including a significant anomalous Hall effect, are observed at temperatures near to T * despite the lack of magnetic ordering 19,41 . The presence of chiral excitonic order would offer an appealing interpretation of the broken time-reversal symmetry and anomalous Hall effect. Second, experiments also report the breaking of threefold rotational symmetry and onset of nematic order around T c ≲ 50 K. Coupling between excitons and CDW naturally results in a phase consisting of nodal d-wave excitons and a nematic CDW, however, our analysis of the free energy found this phase was only ever a local minimum in our model. Coupling to phonons may promote this phase to the dominant ground state, and we leave further examination of this scenario to future work. Third, superconductivity emerges generically as an instability of the TvHS minimal model. However, superconductivity is seen at a much lower temperature scale (T c ≈ 3.5 K) 17,18,20-24 than CDW (T * ≈ 100 K). At these temperatures the Time-reversal symmetry (TRS) breaking and the anomalous Hall conductivity has also been proposed to arise due to a complex CDW state 42,43,52 . Our analysis shows that it is possible for CDW and chiral excitons to coexist, but a key difference between these two states is that chiral excitons break TRS but not translational symmetry, allowing experiment to disentangle the two. To this end, it has been shown that the application of strain and pressure can selectively move the two opposite concavity Fermi surfaces closer or further away from the Fermi level 57 . Moving the Fermi level away from one of the vHS creates a mismatch in the Fermi surface areas, which disfavours the excitonic phase, whereas moving the Fermi level away from the two vHS in a way that keeps the Fermi surface areas roughly equal disfavours charge order, while retaining the tendency to excitonic order. Experimental probes of TRS breaking could be applied in the presence of strain and pressure to disentangle the two phases. Additionally, we suggest that the existence of an exciton condensate should lead to Andreev-like interband tunnelling, an effect which has been used to probe excitonic order in bilayer systems 58 .
Beyond the vanadium metals AV 3 Sb 5 , another possible scenario for topological excitonic condensation is to engineer TvHS in van der Waals heterostructures made from materials with hexagonal symmetry such as graphene or transition metal dichalcogenides (TMDCs) 59,60 . In moiré systems, the reduced bandwidth of the bands near charge neutrality brings the opposite concavity vHS of the valence and conduction bands closer in energy, so that an bias field could feasibly result in a TvHS. The valence and conduction bands can be further aligned through spin splitting the bands via a TMDC layer 61,62 , the effect of which can be tuned via twist angle 63 . Finally, valley polarisation is observed in twisted layered systems [64][65][66] , which could also be exploited to align the valence and conduction bands, e.g., through methods discussed in refs. 67,68 . In the context of layered van der Waals materials, a possible experimental probe would be the enhanced tunnelling between layers induced by excitons, e.g., refs. 69,70 .

Leading instabilities
Our discussion of the leading ordered states follows from the solution of the gap equations for the order parameter vertices We find the mean-field gap equations to be with indices as defined previously: c, d, ± referring to flavour, α to patch, and α denoting the patch connected to α by a nesting vector. To make the equations compact, we have introduced ν = {c, d} with ν = fd,cg. The d-factors are nesting coefficients that characterise the relative strength of the particle-particle and particle-hole divergences, and are defined in the Supplementary Material-we have used notation so that d 0c = 1, d 0d = d 0 , d 2c = d 2 , d 2d = d 3 . The couplings entering the gap equations are understood to inherit scale-dependence from the RG equations for the couplings (5). The eigenvectors for this linear system of gap equations give the possible order parameter structures, and those with the largest eigenvalue are the leading instabilities. The set of Feynman diagrams which generate these flow equations are given in the Supplementary Material.

Landau-Ginzburg analysis
The susceptibility gap equations (10) are insufficient to determine whether order parameters compete or can form a ground state in which multiple orders coexist. Given a set of degenerate or nearly degenerate solutions to the gap equations, we determine which combination of these solutions is the favoured ground state by calculating the Landau-Ginzburg free energy. We employ the mean-field decomposition of the fermions coupled to a combination of order parameter matrices, and integrate out the fermionic degrees of Everywhere we have taken κ = 2; i/rCDW + Φ d+id represents coexistence, and FL the Fermi liquid metallic state. a(i) The iCDW + Φ d+id and b(i) the rCDW + Φ d+id phase diagrams, using the eigenvalues λ Φ , λ CDW as free parameters. In a(ii) and b(ii) the eigenvalues λ Φ , λ CDW are computed explicitly by solution of the RG equations, resulting in a phase diagram as a function of T. The two different phase diagrams result from taking a(i) g 1,c < 0, h 1 < 0, and b(ii) g 1,c < 0, h 1 < 0 with g 3,d = 0; precise values of the initial couplings are given in the Supplementary Material. Article https://doi.org/10.1038/s41467-023-35987-2 freedom, arriving at the free energy Here the full Green's function G À1 ðiω n , qÞ = G À1 0 ðiω n , qÞ + M, ð12Þ comprises the order parameter matrix M = M Φ + M C , The dispersion at each patch is ε 1 ðqÞ = 1 2 q x ðq x + ffiffiffi 3 p q y Þ, ε 2 ðqÞ = 1 4 ðÀq 2 x + 3q 2 y Þ and ε 3 ðqÞ = 1 2 q x ðq x + ffiffiffi 3 p q y Þ. For the two degenerate d-wave excitons, parameterised by Φ a and Φ b , we have Rewriting Tr log G À1 = À F 0 + Tr log 1 where F 0 is the free energy of a free Fermi gas, and using the expansion Tr log 1 + G 0 M À Á = X 1 n = 0 ðÀ1Þ n n Tr ðG 0 MÞ n ð18Þ we evaluate the trace of the first four terms in the expansion, resulting in the free energy stated in the main text. Determining whether the minimum of the free energy can include coexisting C αβν and Φ a , Φ b requires knowledge of the coefficients in this expansion; their calculation is detailed in the Supplementary Material.

Edge states
To demonstrate the presence of edge states in the excitonic phase, we employ a simplified model for numerical diagonalisation, describing a kagome lattice with two-orbital states ν = ± , H = À X hr,r 0 i,ν t ν c y r 0 ,ν c r,ν + X r γ 0 c y r,1 c r,1 + X hr 0 ,ri Δðr 0 , rÞc y r 0 ,ν c r,ν 0 + h.c.

ð19Þ
in which only coupling between nearest neighbours is taken into account. We choose the excitonic pairing function Δðr,r 0 Þ so that the lattice theory possesses an equivalent continuum limit to our field theory description of the three patches surrounding the M points. The gap function is Δðr 0 , r 2 σÞ = 1 ffiffiffi 6 p Δ 0 e iðθ r 0 Àr Àð' + 1Þφ σ Þ , with ℓ = ± 2 equal to the phase winding of the excitonic order around the Fermi surface. The results for ℓ = − 2 are plotted in Fig. 3 in the main text with γ 0 = 2t, Δ 0 = 0.5t, for a ribbon geometry with 60 unit cells.

Data availability
The data produced in this study are available upon reasonable request.